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Abstract 

We study a reaction-diffusion evolution equation perturbed by a space-time Levy noise. The 
associated Kolmogorov operator is the sum of the infinitesimal generator of a Co-semigroup 
of strictly negative type acting in a Hilbert space and a nonlinear term which has at most 
polynomial growth, is non necessarily Lipschitz and is such that the whole system is dissipative. 

The corresponding Ito stochastic equation describes a process on a Hilbert space with dissi- 
pative nonlinear, non globally Lipschitz drift and a Levy noise. Under smoothness assumptions 
on the non-linearity, asymptotics to all orders in a small parameter in front of the noise are 
given, with detailed estimates on the remainders. 

Applications to nonlinear SPDEs with a linear term in the drift given by a Laplacian in 
a bounded domain are included. As a particular case we provide the small noise asymptotic 
expansions for the SPDE equations of FitzHugh Nagumo type in neurobiology with external 
impulsive noise. 
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1 Introduction 

In many areas of investigations, in natural as well as technical and socio-economical sciences, a 
description of phenomena in terms of (partial) differential equations (PDEs) is quite natural and 
has received a lot of attention, also in recent years. However the necessity of taking care of stochastic 
(or random) influences on systems primarly described by (P)DEs in particular through stochastic 
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(P)DEs, S(P)DEs for short, has also came to the forefront of research, see, e.g. |76l H9"] . [3"5].|32j. 
[531 HU [551 H2] In the present paper we concentrate on PDE's perturbed by a space-time noise of the 
additive type. Such SPDE's have been studied extensively particularly in the case of Gaussian noises 
and they have found applications in several areas, from physics to biology and financial mathematics, 
see e.g. [76],[33],[32],[Ml[13[I3[2S[IZl[64]. The extension to the treatment of additive Levy type 
noises (which are more general in the sense that random variables with Levy distributions extend 
Gaussian random variables) is relatively more recent, see e.g. [201 163j . A natural question which 
arises in such extensions from a deterministic description of phenomena to a stochastic one, is in 
which sense one can recover the deterministic description by "switching off' the noise and possibly 
obtain "small noise expansions" around the limit. In the case of SDE's (as opposite to SPDE's) this 
is a rather well studied problem, especially in the case of Gaussian noises and it has also relations 
with the study of the classical limit from quantum mechanics(see e.g. [751 [531 E5J HFJ 030 Q21 (Ml EZJ 
[551 1551 H31 S3])- The infinite dimensional case and the case of SPDEs is less studied, even in the 
case of Gaussian noises, see however [571 HS1 1151 [55J 1SZ1 [55J. Concerning applications, the case of 
stochastic perturbations of the FitzHugh-Nagumo equations of neurobiology and its relations with 
the classical, deterministic FitzHugh Nagumo equations is particularly interesting, due to the fact 
that those equations, and their extensions to the case where the underlying euclidean domain in space 
is replaced by a network, are extensively used in neurobiology, see e.g. J28J [3j EH [76]. In two recent 
papers [5], [J] a systematic study of SPDE's with additive Gaussian noise which includes in particular 
the above stochastic FitzHugh-Nagumo equations, has been given, together with a detailed study 
of the corresponding diffusion expansion around the deterministic limit. One basic difficulty which 
had to be overcome (besides the infinite dimensionality of the stochastic process involved) consisted 
in the non global Lipschitz character of the nonlinear terms. A global Lipschitz condition would 
in fact exclude the interesting case of the FitzHugh-Nagumo equations; similarly other interesting 
equations like those arising in stochastic quantization |12j . hydrodynamics [5J or solid state physics 
( e.g, Allen-Cahn equations), would be excluded. Despite the interest of modeling the noise in such 
systems through a Levy-type one instead of a Gaussian one, which has motivations in all the areas 
which have been mentioned, apparently a corresponding study of asymptotic expansions around the 
underlying deterministic systems has yet to be performed. We shall here adopt the method used in 
[5J to cope with this case. The adaptation involves, in particular, using methods developed by |63j 
to handle stochastic convolutions in the case of Levy noise. Let us note that our results seem to be 
new even in the finite dimensional case, where small stochastic perturbation expansions have also 
not been provided in details for equations of the type we consider. Before we go over to describe 
the contents of the present paper, let us mention that our study of SPDE's with Levy noise can 
also be related to the study of certain pseudo-differential equations with such noises which occur in 
quantum field theory and statistical mechanics (see e.g QJJ1 [H] . Also relations to certain problems 
in the study of statistics of processes described by Levy noises should be mentioned [40l [39] . 
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2 Outline of the paper 



Let us consider the following deterministic nonlinear evolution problem: 

f#(t) = [Attft) + F{4>{t))]dt, t e [0, +00) 

|^(o) = u°, u°eH, [ ' ' 

where A is a linear operator on a separable Hilbert space H which generates a Co-semigroup of strict 
negative type. The term F is a smooth nonlinear, quasi-m-dissipativc mapping from the domain 
D(F) C H (dense in H) with values in H (this means that there exists w G R such that (F — wJ) is 
m-dissipative in the sense of [331 p. 73]), with (at most) polynomial growth at infinity and satisfying 
some further assumptions which will be specified in Hypothesis 13. ll bclow.) Existence and uniqueness 
of solutions for equation (|2.1[) is discussed in Proposition 14.21 below. 

Our aim is to study a stochastic (white noise) perturbation of (|2.ip and to write its (unique) 
solution as an expansion in powers of a parameter e > 0, which controls the strength of the noise, 
as e goes to zero. More precisely, we are concerned with the following stochastic Cauchy problem 
on the Hilbert space H: 



jdu{t) = [Au(t) + F(u(t))]dt + ey/QdL(t), t G [0, +00) 
1m(0)=m°, u° e K , 



(2.2) 



where A and F are as described above, L is a mean square integrable Levy process taking values in 
a Hilbert space U, Q is a positive trace class linear operator from H to H and e > is the parameter 
which determines the magnitude of the stochastic perturbation. The initial datum u° takes values 
in a continuously embedded Banach space K of H . A unique solution of the problem (|2.2j) can be 
shown to exist exploiting as in [24| results on stochastic differential equations (contained, e.g., in 
[32l [33] ) . Our purpose is to show that the solution of the equation (|2.2[) . which will be denoted by 
u = u(t), t G [0, +00), can be written as 

u(t) = 4>(t) + eu\(t) H h e n u n (t) + R„,(t, e) , 

where n depends on the differentiability order of F. The function <f>(t) solves the associated deter- 
ministic problem (|2.1[) . u\(t) is the stochastic process which solves the following linear stochastic 
(non-autonomous) equation 

{dui(t) = [Au 1 (t) + VF(<j>(t))[u 1 {t)]]dt + y/QdL(t), <e[0,+oo) 

\«i(0) = 0, [ ^ ' 

while for each k = 2, . . . , n , Uk(t) solves the following non- homogeneous linear differential equation 
with stochastic coefficients 

(du k (t) - \Au k {t) + VF(0(i))Mt)]] alt + * fc (t)dt, 

\«*(0) = 0. 1 ■ ' 

$fc(i) is a stochastic process which depends on ui(t), . . . , Uk-x(t) and the Frechet derivatives of F 
up to order k, sec Section [5] for details. 
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Let us shortly describe the content of the different sections of the present paper. In Section 3 
we set the basic assumptions needed to perform the construction of solutions and their asymptotic 
expansion. In section 4 we describe the mild solutions to SDE's driven by Levy processes on Hilbert 
spaces, basically following the setting of [63] . Since our expansion will be around solutions of 
the corresponding deterministic equations, we start by presenting results on the latter equations 
(Subsection 4.1). In the Subsection 4.2 we present the setting for the stochastic perturbation, first 
describing the noise. In Section 5 we describe the basic assumptions on the nonlinear term and 
provide its Taylor expansion. Section 6 contains the main results, in particular the construction of 
the expansion, the proof of its asymptotic character and of detailed estimates on remainders, to any 
order. We close with an application to the case of a FitzHugh Nagumo equation on a network. 



3 Assumptions and Basic Estimates 

Before recalling some known results on problems of the types (|2.ip , (|2.2[) , (|2.3[) and (|2.4|) , we begin by 
presenting our notation and assumptions. We are concerned with a real separable Hilbert space, with 
the inner product (■, •). Moreover, in what follows, (B, | ■ |) is a reflexive Banach space continuously 
embedded into if as a dense Borel subset and (K, | ■ | ) is a reflexive Banach space continuously 
embedded in B. On H there are given a linear operator A : D(A) C H — >• H, a nonlinear operator 
F : D(F) C H — >• H with dense domain in H and a bounded linear operator Q on from H to 
H. Moreover, we are given a complete probability space (Q,F, (F t ) t >o, P) which satisfies the usual 
conditions, i.e., the probability space is complete, F contains all P-null subsets of sets in F and the 
filtration (F t )t>o is right continuous. Further, for any trace-class linear operator Q, we will denote 
by TrQ its trace; if / is any mapping on H which is Frechet differentiable up to order n, n G N, 
we will denote by /W, i = 1, ...,n its i-th Frechet derivative and by £)(/W) the corresponding 
domain (for a short survey on Frechet differentiable mappings we refer to Section 4). For any j G N 
and any vector space X, L(X J ;X) denotes the space of j-linear bounded mappings from X 3 into 
X while the space of linear bounded mappings from X into L(X J ;X) is denoted by L 3 (X). We 
denote by | ■ \x the norm on X , by || ■ Hxjpn the norm of any j-linear operator on X and by || ■ \\hs 
the Hilbert-Schmidt norm of any linear operator on X. Finally, for any p > 1, we will denote by 
C?([0, T]; L p (il; X)) the space of X-valued, adapted mean square continuous processes Y on the 
time interval [0, T] such that the following norm is finite 

|||F| = ( sup E|y(^) 1/P <^. 

te[o,T] 

Hypothesis 3.1. 

i. The operator A : D(A) C H — > H generates an analytic semigroup (e tA )t>o, on H of strict 
negative type such that 

\\e tA \\ < e~ ut t > 
\r \\l(h) - e ' - u 

with lu a strictly positive, real constant. 

Moreover, if Ab denotes the part of A in the reflexive Banach space B, that is 
D(Ab) ■= {x e D(A) n B; Ax g B} , A B x = Ax, 
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then A B generates an analytic semigroup (of negative type) e tA]B ,t > on B. 

ii. The mapping F : D(F) C H — > H is continuous, nonlinear, Frechet differentiable up to order 
n for some positive integer n and quasi-m-dissipative, i.e., there exist r\ > such that 

(F(u) ~ F(v) - T)(u — v),u — v) < 0, for all u, v G D(F). 
Hi. If Fjp , j = 1, ... ,n denotes the part of F^ in B, that is 

D{F { B 3) ) := {x e D(F®) fl K;F B 3 \x) G , F^ 5 (x) = 

t/ien i/iere exists a reflexive Banach space K densely and continuously embedded in B which 
makes the following assumptions satisfied: 

(a) there exists a positive real number 7 and a positive natural number n such that: 

\F B (u)\ B < 1 (l + \u\ 2 ? +1 ), u€K, 

(b) for some n and any u £ D(F^),i = l,...,n, there exist positive real constants 7.;, i = 
1, . . . , n such that 

\\F B l) (u)\\ LHB) < 7l (l + \u\ 2 £ +1 ~ l ) , with n as in (ffi), u e K 

iv. The constants uj, 7] satisfy the inequality to — 77 > 0; this implies that the term A + F is m- 
dissipative in the sense of IMj , f33l p. 73]. 

v. The term L is a Levy process (for example in the sense of \22\ \63f ) on some Hilbert space U ; 
moreover we assume that 

f \y\ m v(dy) < 00, 
Ju 

for all m £ IN, where v is the jump intensity measure introduced in section J^.ffp . 

vi. Q is a positive linear bounded operator on H of trace class, that is Tr Q < 00. 

Example 3.2. Let us give an example of a mapping F satisfying the above hypothesis (in view of the 
application to stochastic neuronal models, which we will present in section 6). Let H = L 2 (A) with 
A C R", bounded and open; set B := L 2( - 2n+1 \A), K := L 2 ( 2n+1 ) 2 (A) and let F be a multinomial 
of odd degree In + 1, n <G N, i.e. a mapping of the form F(u) = <?2n+i(u), where g2n+i(u), 
u G H, is a polynomial of degree 2n + 1, that is, <?2n+i(w) = oq + a±u + ■ ■ ■ + a2 n +iu 2n+1 , with 
a t G R, i = 0, . . . ,2n + 1. Then it is easy to prove that D(F) = L 2 ( 2,l+1 )(A) C L 2 (A),n > 0, 
D(F) = L 2 (K) ~ H,n = and (by using the Holder inequality) D(F^) = L 2{ - 2n+1 ~' l \K). Moreover, 
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it turns out that, for any u € K, FW(u) can be identified with the element gm{u) (both in D(F) 
and K). Consequently, 

/ p \ V(2(2n+1)) 

\F(u)\ B ={ I \g 2n+1 (u(0)\ 2{2n+1) dq 

l/(2(2n+X)) 



A 



<c 2n+1 (i + ^K)| 2(2 "+ 1)2 d£ 



C 2n+1 (l + \u\% +1 ) 



and, similarly, 



\y {3) F{u)\ L3{K , B) < C 2n+1 ^(l + M^ +1 - 1 ) 

= C 2 „_i(l + ./ 0.1 '"• 

Hence i 7, satisfies Hvpothesis l3 . f I (ii) . (hi). Further, in the case g3(u) = — it(it— l)(it— £), < £ < 1 the 
corresponding mapping F coincides with the non linear term of the first equation in the FitzHugh- 
Nagumo system (see Example 16.41 below). 

4 Mild solutions to SDE's driven by Levy on Hilbert spaces 

In this section we basically use the setting of [55] . 
4.1 The deterministic case 

Let Ao be a densely defined linear operator on a Banach space B, with domain D(Aq). Let assume 
that the differential equation 

dy A 

= AoV (4.1) 
2/(0) =!,„£ D(A ) 

has a unique solution y(t), t > 0, y(t) € B. The equation being linear, we have 

y(t) = S(t)y , t > 0, 

with S(t) a linear operator from D{Aq) into B. If for each t > 0, has a continuous extension to 
all of B, and for each z € £>, t — > S(t)z is continuous, then one says that the Cauchy problem (|4.ip 
is well posed, t — ► S(t)z, defined then for all z £ B, is called a generalized solution to (|4.1[) . One 
has: that (S'(t)) t>0 is a Co-semigroup: 

i. S(0) = 1, S(t)5(«) = S(t + s), t, s > 0, 

ii. \S(t)z — z\ B — >• as t i 0, for every z e £>. 

Let -D(A) be the definition domain of the generator ^4 of S(t). We have D(A) D -D(^4o) and ^4 is an 
extension of ^4q- Moreover, see, e.g. ([63], Theorem 9.2): 
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i. \S(t)z\ B < e ut M \z\ B , for some w, M > Vz G B, Vi > 0, 

ii. A is closed and for any z G / > one has S(t)z G -D(A) and ^S(t)z = AS(t)z = S(t)Az. 
In particular for z = y G -D(j4), i — ^ S(t)y solves (|4. 1[) with A replacing A . 

Now, let i7 be a Hilbcrt space such that B C H, with a dense, continuous embedding, B being a 
Borel subset of H . Let ^ > be if- valued and continuously differentiable. Then the "variation 
of constants formula,, 

t 

y(t) = S(t)y +Js(t- s)V(s) ds, t>0 (4.2) 
o 

solves 

f t (t)=Ay( t) + m (43) 
2/(0) = y eH. 

In general, whenever the integral in (|4.2p has a meaning for a given yo in H, one says that (|4.2p is 
a mild solution of 

J ^ {4A) 

I 2/(0) = 2/o- 

The formal definition of mild solution is given below. 

Definition 4.1. Let yo G iT; we say £/iai £/ie function 4> : [0, oo) — > if is a mild solution of equation 
(|2.ip i/ is continuous (in t), with values in if and it satisfies: 

4>{t) = e tA y + f e('- s ^(s) ds, t G [0, +oo), (4.5) 



with the integral existing in the sense of Bochner integrals on Hilbert spaces. 

In the case of ijj being substituted by a mapping F : D(F) C if — >• if satisfying the assumptions 
given in Hypothesis 13 . 1 1 we have the following result. 

Proposition 4.2. Under Hypothesis 13.11 there exists a unique mild solution <f> = <f>(t),t G [0, oo) of 
the deterministic problem 

2/(0) = j/o G D(A ) 
such that 

W)\h <e- 2 ^-^\u Q \ H , t>0. (4.7) 
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Proof. The proof of existence and uniqueness can be found, e.g in [321 Theorem 7.13, p. 203], 
while estimate (|4.7[) is a direct consequence of the application of Gronwall's lemma to the following 
inequality 



Remark 4.3. It can be shown that, under Hypothesis 13.11 there exists a K -continuous version of 
the unique solution of equation (|4.5[) such that, for any T > 0, p > 1 



(see f3Sl Section 5.5.2, Proposition 5.5.6]). Hence, in the following, by <p we will always understand 
this K -valued version of the solution of (|2.1[) . 

4.2 The stochastically perturbed case 

Let G be a linear operator from a Hilbert space U into a Hilbert space H . Let S(t) be a Co-semigroup 
on the Hilbert space H. Assume the generator (A, D(A)) of S(t) in H is almost m-dissipative (i.e. 
[(Al — A) + rj\H = H for any A > and some r\ E R: [|63j. p. 180]; this is equivalent to quasi 
m-dissipative in the sense of [33J P- 73].) Assume B C H as in (|4. 1 [) and that the restriction Ab of 
A to B is also almost m-dissipative. Let L be a square-integrable mean zero Levy process taking 
values in a Hilbert space K. I.e. L = (L(i)J >Q takes values in a Hilbert space if, has independent, 
stationary, increments, one has L(0) = 0, and L(t) is stochastically continuous (see [S3], Def. 4.1, 
p. 38). Let Q be the covariance of L. Then Q^(K) is the reproducing kernel Hilbert space (RKHS) 
of L, assume Q*(K) is embedded into U. 
We recall the following basic notions and results. 

Definition 4.4. Let v be a finite measure on a Hilbert space U such that ^({0}) = 0. A compound 
Poisson process with Levy measure (also called jump intensity measure) fi is a cadlag Levy process 
L satisfying 



^-Mt)\ 2 H = 2{A4>(t),4>(t))dt + 2{F{<p{t)), 4>{t)) 
at 

< -2(u-r,)\<j>(t)\ 2 H- 



□ 



sup \<f>(t)\ 



< 00. 



te[o,T] 





S(/7) being the a— algebra of Borel subsets ofU. 



Given a Borel set I separated from 0, write 



7r / (*) = X;X/(Ai(a)) > 



t > 0. 



s<t 



8 



The cadlag property of L implies that 717 is Z + -valued. We notice that it is a Levy process with 
jumps of size 1 and thus, a Poisson process (see [[B3J, Proposition 4.9 (iv)] for more details.) We 
also have that E7T/(i) = tE7T/(t) = tv(I), where v is a measure that is finite on sets separated from 
0. We shall write 

Lj(t) = 5>j(AL(s))AL(s). 

Then Lj is a well-defined Levy process. The theorem below provides the corresponding Levy- 
Khinchine decomposition: 

Theorem 4.5. i. If v is a jump intensity measure corresponding to a Levy process then 

{\V\1 A l)i/(dy) < 00. 

r 

ii. Every Levy process has the following representation: 

L{t) :=at+^QW(t)+jr(L Ik (t)-t f yu(dy)) + L Io (t), 
k=i ^ Jlk ' 

where Iq := {x : \x\jj > ?'o}j Ik '■= {% ■ ru < \x\xj < Hc-i}, (ffc) is an arbitrary sequence de- 
creasing to 0, W is a Wiener process, all members of the representation are independent 
processes and the series converges P — a.s., uniformly on each bounded subinterval [0,oo). 

In the following (see Hypothesis 13. ip . with no loss of generality, we assume that 

/ 1/Kdtf) = 0. (4.8) 

We also assume throughout that the Levy process is a pure jump process, i.e. a = and Q = and 
that 



\y\ m v(dy) < 00, for all m € N, (4.9) 
which leads to the representation 

00 

L(t) = J2 L i*(t)+Li (t), 

k=l 

in view of assumptions (|4.8p and (|4.9[) . 

t 

Let -La(£) = / S{t — s)v / Qdi(s), t > 0, be the Levy Ornstein-Uhlenbeck process associated with 


S, y/Q, L, assumed to exist and have a cadlag version in B (the latter is satisfied if B is a Hilbcrt 
space K and Sit) is a contraction on K), see e.g. ([63], p. 155), or 5 is analytic and L takes values 
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in £>((— A) a ) for some a > 0; see, e.g. ([S3], P- 155). Assume F is an operator on H (possibly 
nonlinear, nor everywhere defined) satysfying Hypothesis 13.11 

An adapted B- valued process X is said to be a cadlag mild solution to 

f dX(t) = AX(t) dt + F(X(t)) dt+^Q dL(t) 

\ X(0) = x£ D(F) 1 ' ' 

if it is cadlag in B and satisfies, P-a.s., the equation 
t 

X(t) = S(t)x + J S(t-s)F(X(s))ds + L A (t), t > 0, with X(s) £ D(F) for s > ([53], p. 182). The 
o 

formal definition of mild solution for the stochastic problem (|4.10[) is given below; next we recall the 
definition of stochastic convolution and we list some of its properties. 

Definition 4.6. Let u° £ K. A predictable H -valued process u := (u(t))t>o is called a mild solution 
to the Cauchy problem (|2.2|) with initial condition u° £ D(F) if for arbitrary t > we have 

otA.,0, f t ,{t-a)A J? (,.( B \\ da , _ f Jt-s)A 



u(t) = e tA u { >+ e {t - s > A F(u(s))ds + e e ( ' yjQdL(s), P-a.s. 
Jo Jo 

La{J,) '■= Jq e <it ~ s ' >A \/QdL(s) is called a stochastic convolution and under our hypothesis it is a well 
defined mean square continuous S^t-adapted process with values in B and cadlag trajectories (see 
e.g., JEM, Proposition 9.28, p. 163). 

The first integral on the right hand side is defined pathwise in the Bochner sense, P-almost surely 
For further use, in the following we introduce some additional condition on the stochastic con- 
volution: 

Hypothesis 4.7. The stochastic convolution L-A{t),t > introduced in Definition 14.61 admits a 
K-valued version such that, for any T > 0, it satisfies the following estimate 



E sup \L A (t)\% <C T (4.11) 
\te[o,T] J 

for every m € hi and some positive constant Ct (possibly depending on T) . 

Example 4.8. Let us give an example for the setting (H, B, K, L, A,Q) where La is well-defined 
and Hypothesis 14.71 is satisfied. This example is related to the application to the stochastic FitzHugh- 
Nagumo model which we discuss in Example 16.41 Let H, B, K be as in Example 13.21 Let A = A be 
the Laplacian in L 2 (K) with Neumann boundary conditions on the boundary d A of the bounded 
open subset A of R". Let Q be a bounded trace class operator commuting with A and L be a Levy 
process such that the corresponding measure v satisfies 

M^,3 2(2„ + i)f(da;) < oo for all m £ N, 

L*(A) 

where iy /3 ' 2 ( 2n+1 - ) is a fractional Sobolev space with given (3 > 0. Finally, let (Ai$, D(Ai$)) denote 
the generator of the heat semigroup with Neumann boundary conditions on L 18 (A). By |32l Ap- 
pendix ] La(£) € D((— ^2(2n+i))' 7 > 0; in particular Lj^{t) £ K, La being in addition a Levy 
process. This implies the bound in Hypothesis 14.71 
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The next result concerns the existence and uniqueness of solutions for the stochastically perturbed 
problem. Moreover, we shall use Hypothesis 14.71 above concerning the Ornstein-Uhlenbeck process 
associated with e tA , ^/Q and L in order to prove a useful estimate on the solution. 

Theorem 1. Assume that A and F satisfy Hvvothesis YS.ll Assume that A and Q satisfy Hypothesis 
14.71 Then there exists a unique cadlag mild solution of for any x G B. For each x G H 

there exists a unique generalized solution for 10\ l (in the sense that 3(X n ) ne ^, X n G B, unique 
cadlag mild solutions of (|4.10[) with X n (0) = x s.t. \X n (t) — X(t)\ H — > uniformly on each bounded 
interval). Moreover defines Feller families on B and on H (in the sense that the Markov 

semigroup P(t) associated with X(t) maps for any t > 0, Cb(H) into Cb(H) and Cb{B) into Cb(B) ). 
Moreover, the solution X to J^. Hfy belongs to the space L P (Q; C([0, T]; H)), i.e., is such that 

E( sup \X{t)\ p H \ <+oo, (4.12) 

\t£[0,T] J 

for any p G [2, oo). 

Proof. The first part of the result is proven in ([63], Theorem 10.14). We only have to prove the 
estimate ()4.12|) . Let z(t) := X(t) — L A (t); then it is not difficult to show that z(t) is the unique 
solution of the following deterministic equation: 

(z'(t)=Az(t)+F(z(t) + L A (t)) 
\z(0) = u° 

with z'(t) := f t z{t). 

With no loss of generality (because of inclusion results for L p -spaces with respect to bounded 
measures) we can assume that p = 2a, a G IN. Now combining condition (i) with (i) in Hypothesis 
(|3.ip and recalling Newton's binomial formula we have: 



-\z(t)\% a = 2a(z\t),z(t))\z(t)\%- 2 =2a(Az(t) + 

< -2aoj\z(t)\% + 2a(F(z(t) + L A (t)), z(t))\z(t)\ 2 "- 2 

< -2a(u> - v)\z(t)\ 2 H a + 2a\F(L A (t))\ H \z(t)\ 2 H a - 1 

< -2a(uj - vMt)\ 2 H + 2a^\F(L A (t))\H + C a 2a£\z(t)\ 2 H a , 



(4.13) 



for some constant C a > and a sufficiently small £ > such that —2a(u> — r/) + 2a£C a < 0. Applying 
the previous inequality and Gronwall's lemma we get: 

\z(t)\jf < e (-2a(— I? )+ 5 C a 2a) t | u 0|2a + ^ £ -2a( U - r,) (t-s) | {g) j 1 2a ^ 

S Jo 

Then there exists a positive constant C such that: 

\X(t)\j? < C ^ e (-2a(^)+?C o 2a)t| u 0|2a + 2fl ^ e ~M^-n) (*-») ( s ) ) | ja ^ + |£ A ( t )||« 



(4.14) 
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Since by condition (iii) in Hypothesis 13. 11 the restriction of F to K has (at most) polynomial growth 
at infinity in the if-norm and, by the assumption on L A (t) made in Hypothesis 14.71 L A takes value 
in K, for any a £ hi we have: 

\F(L A (t))\jf < C a>m (l + \L A (t)\^) 2a < C a , m (l + \L A (t)\j? m ), 

for some positive constant C a _ m depending on to and a. Moreover, we observe that, again by 
Hypothesis 14. 7\ it holds that 

e( sup \L A (t)\% m ) <c; m>r , 

where C' a m T is again a positive constant depending on to, a and T; hence 



sup / e 

te[o,T] Jo 



-2a(cu — 7]){t— s) 



\F(L A (s))\ 2 H a ds 







< CE 


sup / 




te[o,T] Jo 



< CE 



■t rt 

sup / e- 2a( "-" )( *- s) ds + C;. m / e^^ds 

t£[0,T]Jo ' Jo 



<C, (4.15) 



for some positive constants C, C depending on a, m and T. Consequently, putting together inequal- 
ities (|4.14j) . (|4.15|) . we obtain 



e sup \x{t)\% <c\ u y H a + c 

\te[o,T] I 



for some positive constant C, so that the proposition follows. 



□ 



5 Properties of the non-linear term F and Taylor expansions 

In this section we study the non-linear term F in order to write its Taylor expansion around the 
solution cj)(t) of (|4.5[) with respect to an increment given in terms of powers of e. In order to do that 
we recall some basic properties of Frechet differentiable functions. 

Let U and V be two real Banach spaces. For a mapping F : U — > V the Gateaux differential at 
u G U in the direction h € U is defined as 

s->0 s 

whenever the limit exists in the topology of V (see for example [52l p. 12]). 

We notice that if \7F(u)[h) exists in a neighborhood of uq £ U and is continuous in u at uq 
and also continuous in h at h = 0, then WF(u)[h] is linear in h (see for instance [52l Problem 1.61, 
p 15]). If VF(mq)[/i] has this property for all uq € Uq C J7 and all ft- G C/ we shall say that F 
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belongs to the space G 1 (Uq;V). If F is continuous from U to V and F G G 1 (Uq; V) and one has 
F(u + h) = F(u) + VF(u)[h] + R(u, h), for any u G Uq with: 



im 1 \ u ' lv = 5.1 
\h\ v 

with | • |y and | ■ |j/ denoting respectively the norm in 1/ and U, then the map /i — ¥ X7F(u)[h] is a 
bounded linear operator from Uq to V, and V.F(u)[/i] is, by definition, the unique Frechet differential 
of F at u G Uq with increment ft- G U. The function h) is called the remainder of this Frechet 
differential, while the operator sending h into VF(u)[h] is then called the Frechet derivative of F at 
u and is usually denoted by F'(u) (see for instance [52l pp. 15-16, Problem 1.6.2 and Lemma 1.6.3]). 
We have then VF(u)[h] = F'(u) ■ h, with the symbol • denoting the action of the linear bounded 
operator F'(u) on h. 

The mapping F'(u) is also called the gradient of F at u (see for example [52j p. 15]) and it 
coincides with the Gateaux derivative of F at u. We shall denote by (Uq,V) the subset of 
G 1 (Uq,V) such that the Frechet derivative exists at any point of Uq. Similarly we introduce the 
Frechet derivative F"(u) of F' at u G U. This is a bounded linear map from a subset D(F') of 
U into L(U, V) (L(U, V) being the space of bounded linear operators from U to V). One has thus 
F" G L(U, L(U, V)). If we choose h,k G U then F"(u) ■ k G L(U, V) and (F"(u) ■ k) ■ heV. The 
latter is also written F"(u) h k or F"(u)[h, k]. The mapping F"(u)[h, k] is bilinear in h, k, for any 
given u G D{F") and it can be identified with the Gateaux differential F{u)[h, k] of \7F(u){h] 
in the direction k, the latter looked upon as a map from U to L(£7, V). Similarly one defines the 
j-th Frechet derivative F^\u) and the j-th Gateaux derivative VF^'(u)[/ii, . . . , /ij]. The function 
FV>(u) acts j-linearly on h\,...,hj with ft,, G U for any i = 1, . . . , j. Let Uq be an open subset of 
U and consider the space 5W{U ,V) of maps F from {/ to V such that F^) (u) exists at all u G Uq 
and is uniformly continuous on Uq. The following Taylor formula holds for any u, h G U for which 
F(h) and -F(it + /i) are well defined (i.e. h and u + h are elements of D(F)), and j = 1, . . . , n + 1 
with u G n? =1 3V)(U ,V): 

F{u + h)= F(u) + S7F(u)[h] + iv (2) F(u)[/i, h] + ■ ■ ■ + ^V (n) F(u) [h, . . . , h] +R ( - n \u; h) , (5.2) 

where \R^(u;h)\ u < C u>n -\h\^ for some constant C Ui „ depending only on u and n (see for example 
[52 Theorem X.1.2]). 

Now let us consider the case U = H 7 with H being the same Hilbcrt space appearing in problem 
(|2.ip . Let F be as in Hypothesis 13.11 and set Uq — D(F). Let us define for < e < 1 the function 
h(t), t > 0: 



h(t) = J2e k Mt) + r {n) (t;e), 



k=l 



where the functions Uk(t),k = l,...,n and r^ n \t;e) are p-mean integrable continuous stochastic 
processes with values in H, defined on the whole interval [0, T] for p G [2, oo). Moreover we suppose 
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r («)(. ;e ) = (s n ), i.e., 



lim E 



sup 

*e[o,T] 



|rW(i;e)|P 



0, for any T > 0. 



Let be a p-mean integrable continuous stochastic process with values in the Banach space K. 
Then using the above Taylor formula we have 



F(m + h(t)) = F(<f>(t)) + VF(#))[Mi)] + \y {2) F[h(t), h(t)} +■■■ 

+iV(")f(«) [h(t), . . . , &(*)] +i? (n) (<X0; M*)) , 



(5.3) 



and, recalling that for any j = 1, . . . , n, V^F(</>(t)) is multilinear, we have 

W»F(mw),...,h(t)} = 



j-terms 



1 "J 

£ £fel + ... + fc . V O) F(0W)[ufei (t)) _ ^ (t)] + 0j{£ u 3) 



(5.4) 



fciH hkj—j 



where Oj (e nj ) is the contribution to the right member of the above equality coming from the term 
r^ n ' (t; e) and satisfies the estimate 



lim E 

e-s-0 



sup 

te[o,r] 



0, /or any T > 0. 



We notice that any derivative appearing in the member on the right hand side of (|5.4[) is multiplied 
by the parameter e raised to a power between j and nj. 

Taking into account the above equality we can rewrite (|5.3[) as 



F(cP(t) + h(t)) = F(cf>(t)) + J2 £ fc VF(0(t))K(t)] 



E : 

h+]2=2 



k=l 
TJi+h 



2! 



-V^F{4>{t))[ Ujl {t),u h {t)] 



(5.5) 



fcl 



+ E ■ 

iiH hjfc=fc 

+ vWf^WjWt), . . . , Ul (t)] + i?i n) (^); ft(t),e) 



where the quantity Rj n ((f>(t);h(t),e) is given in terms of the derivatives of i 7, with the parameter 
e raised to powers greater than n, in terms of the n-th remainder R^ n \(f>(t);h(t)) in the Taylor 
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expansion of the map F (as stated in equation (|5.2I) ) and in terms of the remainders Oj{e n3 ), 
j = 2, . . . , n introduced in (|5.4p . Namely, we have: 



n nj 



j— 2 iiH Hj— n+1 

n 

3=2 



(5.6) 



R( n >((f>(t); h(t)) being as in (|5.2j) (with u replaced by <^>). In this way equation (|5.5p can be rearranged 
as 

F(<t>(t) + h{t)) 



F(o(*)) + E £i E -v (j) ^(t))KW,...,^.(i)] 

j=2 \i!H hi i= i ■ 7 " 

^ n) (0(t);ft(t), £ ). 



(5.7) 



Lemma 5.1. Let R^ 1 be as in formula (|5 .6[) . Then for all p £ [2, oo) tm<i T > £/iere exists a 
constant C > 0, depending on \4>\k, ■ ■ ■ , \u n \ii, V^i* 1 , . . . ,V( n 'F,p,n, such that: 



E 



sup |4 n) («^);Mt) l£ )rH 

te[o,T] 



/or oK < e < 1. 

Proof. First of all we notice that 

meaning that 



X>( £ ^) = o( £ 2 »), 

3=2 



E°( £,y ) 

3=2 



< C„e 2 ™, e -> 0, 



(5.8) 



for some constant C„ > 0. Now since: 

n nj -, 

R[ n) (4>(t); h(t), e) = E E £il+ '" +t3 -V W) F(0(t))[ U4l (t), . . . , Uij (t)] 



j— 2 ii+...+ij=Ti+l 
n 

^ 0j ( e ^) +J R(")(^);Mt)), 

3=2 
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using the estimate given in condition (3.b) in Hypothesis 13. II and (|5.8[) . for e £ (0, 1] we have 
\R^(mh(t),s)\ p H 



maxJ|V«F(^(t))|| L3 - w J I 



pin 



I'' 
I if 



(5.9) 



< CWe(" +1 ^.maX n [ 7 f (1 + W^Y] (E l«i(*)ISr) 

+ a i e 2 " p + C( 2 )|i?^(^);/.(t))|? f 

<C ne ("+ 1 )" + cW|ieW(^(t);/ l (t))ISr. 

where p -,Cn]>,C^l are constants depending only on n,p and the constant C n in (|5.8p while 

C„ is a suitable positive constant depending on p, n, maxj = i „ 7^(1 + |0(^)lx ~ 3 ) p ili being the 

constants appearing in Hypothesis 13. 1[ condition (3)) and \ui(t)\ p H , i = 1, . . . ,n. We notice that the 
above inequality follows by recalling that the deterministic function (j){t) is bounded (in the TJ-norm) 
(see Proposition ^. 2[) . 

Now by the bound on in the equation (|5 .2[) we have that 



(n+l)p 



\RW(<f>(t);h(t))? H <C n \h(t)\% 

with C n depending on </>(i) and n but independent of h(t). Since h(t) = Y^k=i £ku k{t) + r^ n '{t]e) 
with \A n \t;e)\ < C n e n+1 for some C n , then: 



\R^m); h(t))\ p H < e^ +1 ^C n , p (\ Ul (t)\ H , \u n (t)\ H ) 

with C niP = C n , P (\ui(t) \ H ,.. . , \u n (t)\ H ) independent of e. 
Hence by (|5.9|) and (|5.10|) we have that 



(5.10) 



sup |4 n) (0(<);MO,£)l?/ 

*e[o,T] 



< C' n e n+1 , 



where C' n := C' n {p, V^F, . . . , V^F, |</>|#, . . . , \u n \n ) is independent of e. This gives the lemma, 
with C = C' n . 

□ 

As we said before, we want to expand the solution of the equation (|2.2|) around </>(£), that is we 
want to write u(t) as: 

u(t) = <j){t) + e Ul (t) + ■ ■ ■ + e n u n {t) + R n it, e), (5.11) 



1G 



(with the term Rn(t, e) = 0(e n+1 )), for any t > 0), where the processes (v,i(t))t>o, i = 1, • • • , n can 
be found by using the Taylor expansion of F around <p(t) and matching terms in the equation (|2.2p 
for u. Given predictable TJ-valued stochastic processes w(t), V\(t), . . . ,v n (t) let us use the notation: 

k 

S*Mt)) [«!(*),..., E vWfWt))K(t),...,,,i, (5.12) 

3=2 ii H Mj=fc 

with ix, . running from to fc and the given restriction i\ + ■ • • + i n = k. With the above 

notation the processes u\(t), . . . ,u n (t) occurring in (|5.11[) satisfy the following equations: 

Um(t) = [A Ul (t) +VF{(j)(t))[u 1 (t)]}dt + y/QdL(t), 
\u 1 (0) = 0, 

and 

idu k {t) = [Au k {t) + VF(cf>(t))[u k (t)}]dt + $ k (t)dt, 
\u k (0) = 0, 

with 

* fc (t) := # fc (0(t)) [ui(t), • • ■ , u fc _i(t)] , fc g N, n > fc > 2 . (5.14) 

Notice that while Ui (t) is the solution of a linear stochastic differential equation (with time dependent 
drift operator A + WF((p(t))), the processes U2, ■ ■ ■ , u n are solutions of non-homogenous differential 
equations with random coefficients whose meaning is given below. 

Definition 5.2. Let 2 < k < n. Then a predictable H-valued stochastic process Uk = Uk(t) ,t > is 
a solution of the problem (|2.4[) (i.e. (|5.13[) ) if almost surely it satisfies the following integral equation 

«fc(*) = / e^ A VF((f>{s))[u k {s)}d S + f $ fc (s)ds, t > , 2 < k < n, 



with <f> as in Provosition 14.21 and $fc as in (|5.12p and (|5.14|) . 

In the following result we estimate the norm of $fc in by means of the norms of the Gateaux 
derivatives of F and the norms of Vj(t), j = 1, ...,k — 1, where fj(i) are i?- valued stochastic 
processes. 

Lemma 5.3. Lei us fix 2 < k < n; let w(t) and Vi(t), . . . ,v k —i(t) be respectively a K-valued 
process and H-valued stochastic processes. Then ^ k {w{t)) [vi(t), . . . , Vk-i(t)] as in (|5. 12|) satisfies 
the following inequality 

|**M*)) M*), ■ • • , < chW|^fc 2 (fc + Mi)^ 1 + ■ • • + l^-xWIir 1 ), 

where C is some positive constants depending on k and the constant jj, j = 2, . . . , k introduced in 
Hypothesis 13.11 
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Proof. We have 

\* k (w(t))[v 1 (t),...,V h - 1 (t)]\H 



E E 

j=2 iiH M 3 =fc 



^F(w(t))[v il (t),...,v i .(t)] 



^E E 

3=2 iiH \-ij = k 



V«F( W (i))K (<),..., «,,(*)] 



(5.15) 



and using the assumption (3) in Hypothcsis l3.il we get 



k j 

\Mt)\H<fl E ^\\VF^(w(t))\\L H H } i[K(t)\H 
j=2ii+...+i j =k J ' 1=1 

<E \i^ + \^)wr' 3 E EkwIh 

j'=2 iiH hij =& ( = 1 

fc / fc-1 N 

< E + M*)i*r -i E u+E 1 

i=2 h+—+ij=k V J=l / 

fc / fc-1 \ 

< E + M*)i*r~ J '* 2 U+E I^WIir" 1 



(5.16) 



3=2 



1 = 1 



fc-1 



< 



c(i + \w(t)\™- 2 )k 2 fc + ^hWIff 1 > 



i=l 



for some positive constant C, from which the assertion in Lemma 15.31 follows 



□ 



Remark 5.4. Notice that by Lemma f5.3[ if vi, . . . , Wfc-i are p-mean (p € [2,oo)), integrable con- 
tinuous stochastic processes then the same holds for 

6 Main results 

Proposition 6.1. Under Hypothesis 13.11 the following stochastic differential equation: 

(6.1) 



d Ul (t) = [Aui(t)+VF(4>(t))[ui(t)]]db + y/QdL(t), te [0,+oo) 
«i(0) = 0, 



has, with <p as in Proposition 14.21 a unique mild solution satisfying, for any p > 2, the following 
estimate: 



E 



sup \ui(t)\ p H 
te[o,T] 



< +oo, for any T > 0. 



(6.2) 
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Proof. First we show the uniqueness. Let us suppose that w\(t) and w% (t) are two solutions of (|6.1I) . 
Then by Ito's formula we have: 



d\ Wl (t) - w 2 (t)\ 2 H = (A( Wl (t) - w 2 (t)),w 1 (t) - w 2 (t)) dt 

+ (WF(cb(t))[ Wl (t) - w 2 (t)] , Wl (t) - w 2 (t)) dt, 

so that, by the dissipativity condition on A and the estimate on VF in Hypothesis 13. 11 (3), we have 

d\wx(t) - w 2 (t)\ 2 H < -Lu\ Wl (t) - w 2 (t)\ 2 H + 7 i(l + |0|£ _1 )Mt) - w 2 (t)\ 2 H . 

Now uniqueness follows by applying Gronwall's lemma. 

As far as the existence is concerned, we proceed by a fixed point argument. We introduce the 
mapping F from £ P {Q- C*([0, T];H)) into itself defined by 

T(w(t)) := f e^ s ~> A VF^( S )))[w(s)}ds + L A (t). 



We are going to prove that there exists T > such that T is a contraction on L' p (il; C([0, T]; H)). 
In fact, for any v, w G L p {Vt; C([0, T]; H)) we have, for any < t < f : 



|r(«(t))-rMt))|p 



sup 

te[o,T] 



=(*-«)^VF(0(a)))[t;(«) - w(s)]ds 



h 



< E 



< E 



< E 



sup 

te[o,f] Jo 



e^ A \\l {H) \VF(0( S )Ms)-w(s)}\ p H ds 



sup \VF(<j>(s))[v(s)-w(s)]\ p H 
se[o,f] 



\\e^ A \\% H) ds' 



sup |u(s) -w(s)|# 

s6[0,f] 



7f(i + i^)irT^(i-^ pf ) 



1 - e-"v T 



where we used condition (Hi) in Hypothesis 13.11 for the third inequality and Proposition 14.21 for the 
last inequality. Then if T is sufficiently small (depending on tu,p, 71, </>), we see that T is a contraction 
on &(Sl;C([0,f\;H)). 

By considering the map T on intervals [0, f ], [f, 2T] , . . . , [(AT - 1)T, T], f = T/N, N G N, we 
have that T is a contraction on £ p (f2; C([0, T]; iJ)) and hence we have the existence and uniqueness 
of the solution for the equation (|6.1[) in the space L p (il; C([0, T]; H)) for any p G [2, 00). 

Let us now consider the estimate (|6.2j) . We write the ltd formula for the function | • \ jf H, applied 
to the process X. To this end, we recall the expressions for the first and second derivatives of the 
function H(x) := |x| 2a ,a G N. 
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We have 

VF(a:) = 2a\x\ 2( - a ~^x 

iTr(QVF 2 (x)) - aTr(Q)| a; | 2 ( Q - 1 ) + (a - l)a\x\ 2 ^\y/Qx\ 2 . 

Moreover (see, [53] and [3(3), we recall that ltd formula implies: 

dF(u(t)) = WF(u(t-))du{t) + ^Tr(QVF 2 (u(t-)))du + d[u](t) 

Although our computations are only formal, they can be justified using an approximation argument. 
By condition (iii) in Hypothesis 13.11 we have for all points in the probability space and p = 2a with 
a G N: 

dMt)| 2a ^ 2a( Ul (t-),d Ul (t)) H \ Ul (t)\ 2a ~ 2 + aTT(Q)\ Ul (t^)\ 2( - a -^dt 

+ (a- l)aMH| 2 ( a - 2 )|yQ Wl (t)| 2 dt + d[«i(t)](t). (6.3) 
By the dissipativity of A + F, the first term in the above inequality is estimated by 

( Ul (t-),d Wl (t)) ff K(t)| 2a - 2 
= {Amitlmmmit)^- 2 + (VF(m)[Mt)},Mt))\Mt)\H~ 2 + (VQdiKt-W^Mt)!*?- 8 

< - w|m(t)||? + 2a 7 (l + lu ^" 1 )!^^)!^ + 2a(^/QdL(t), ui^M*)^" 1 

(6.4) 

where uj := uj — 7(1 + \u 

Moreover, the second and third term in (|6.3[) can be estimated in the following way: 

aTr(Q)\ Ul (t)\ 2 ^ + (a - l)a\ Ul (t-)\ 2 ^\ ^Q Ul (t)\ 2 < C a (eTv 2a (Q) + ^i(i)| 2 °), (6.5) 

for any e > 0, where we used the elementary inequality a6 2 ^ a_1 - ) < C a (ea 2a + \b 2a ), with C a being 
a suitable positive constant. Therefore 

|ui(t)|l? < -2a {^~^f) ^Vi(s)|^d s + 2a£( v ^dL( s ), Wl ( s ))| Ul ( s )| 2 ^- 1 

+ C a eTr(Q) 2a T + f TrQd|i|(s) 
Jo 

and 



C 

Esup|ui(t)|^ a < -(2auT -)Esup \u x {t)\ 2 £ 

t<T e t<T 



2a E sup 

t<T 



t 

Qdi(s) ) ui(s))|wi(s)| 2 |- 1 ds 



CatT + T [ TrQ|a;| 2 ^(da;), 
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where we used the relation 



t<T 



Esup[ui](t) < E / Tr(Q)d[L](t) = E / Tr(Q)d(L)(t) = T Tr(Q)\x\ 2 v{dx). (6.6) 



H 



By the Burkholdcr-Davis-Gundy inequality, see.e.g, ([53], p. 37, HU05]) applied to 



M(t) := / (V QdL(s),ux(s))\ui(s)\x 
Jo 



there exists a constant c\ such that 

Esup / 
t<T Jo 



QdL(s), Ul (s-))\ Ul (s)\%- 



<ciE 



QdL(8),ui{a-))\ui(a)\$- 



(T) 



1/2 



< ci E (sup \ Ul (t)\ 2a [ Tr(Q)d[L](s)\ 

\t<T Jo J 

< deEsup|«i(t)|H / TrQ\x\ 2 v(dx), 



t<T 
,2 j_ f-, Ia c \u2 



where we used the elementary inequality ab < ear + (l/4e)6 , e > 0. Collecting the above estimates 
we obtain 



Esup \ux(t)\% < -{2a£j - — )TEsup \ Ul (t)\j? 



t<T 



t<T 



+ 2 Cl eEsup|ui(i)|^ + (%- + i)t [ TrQ\x\ 2 [t(dx) + C a eT, 



Hence 



sup \ Ul (t)\j? 

te[o,T] 



<T r" -2a(cD-c a /e)T , ^ 



where C a ,T is a positive constant and (|6.2[) follows. 



□ 



Theorem 6.2. Lei us /ix 2 < k < n, assume that Hypothesis \3.1\ holds, and let U\ be the solution 
of the problem (|2.3j) . Suppose moreover that Uj is the unique mild solution of the following Abstract 
Cauchy Problem (ACP): 



duj(t) = [Auj(t) + VF(<j>(t))[uj(t)]]dt,+®j(t)dt 
ttj(0) = 



for j = 2, . . . , k — 1 satisfying: 



sup K'(*)|?j 
te[o,T] 



< +oo, T > 0, for any p <G [2, oo); 



(ACPj) 



(6.7) 
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then there exists a unique mild solution Uk(t) of the following non-homogeneous linear differential 
equation with stochastic coefficients (in the sense of Definition[ 



(ACP k 



\du k (t) = [Au k (t) +WF((f>(t))[u k (t)}}dt + Q k (t)dt, t £ [0,+oc) 
|u fc (0) = 

and it satisfies the following estimate, for any T > : 



sup \u k (t)\ p H 
te[o,T] 



< +00. 



Proof. We proceed by a fixed point argument, where the contraction is given by 



(6.8) 



e (- t -^ A VF(<j>(t))[y(t)]ds+ / e^ A ^ k (s)d. 



r(»(t)) 



on L p (fl;C([0,T];H)). In fact, arguing as in Proposition 16. II we see that for T £ [0, T] sufficiently 
small, T is a contraction on L p (il; C([0, T]; H)), p £ [2, 00), so that the existence and the uniqueness 
of the solution for QACPkD follows. 

Let us consider the estimate (|6.8p . By the condition (iv) in Hypothesis 13.11 we have, for p = 2a 
with a £ N (and all points in the probability space) : 



-\u k (t)\l a = 2a(Au k (t),u k (t))\u k (t)\l a - 2 + 
+ 2a{<$> k {t),u k {t))\u k {t)\ 2 "- 2 

< -2aoj\u k (t)\j? + 2a 7 (l + \u°\ K )\u k (t)\% + 2a\<f k (t)\ H \u k (t)\^' 1 

< -2aCj\u k (t)\j? + C a \$ k (t)\j?, 



(6.9) 



where uj :=a; — 7(1 + |u°|if) as in the proof of Proposition (|6.1I) . By the assumption (|6.7|) made on 
Uj(t),j = 1, . . ., k — 1 and Lemma 15.31 we have that: 



SUP \$ k (t)\% a 

te[o,T] 



<C' a , T > 0, 



so that taking the expectation of inequality (|6.9[) and applying Gronwall's lemma (similarly as in 
the proof of Proposition 16. ip we obtain: 



sup \u k {t)\j? 
te[o,T] 



< C'e- 2aGiT < C a 



where C a is a positive constant, and the theorem follows. 



□ 



We are now able to state the main result of this section: 
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Theorem 6.3. Under Hvvothesis YS.ll the mild solution u(t) of (|2.2[) (in the sense of Definition 
can be expanded in powers of e > in the following form 

u(t) = (f>(t) + EUxit) H \-e n u n (t) + R n (t,e), n e N, 

where u\ is the solution of 



dui(t) = [Aux{t) + VF(0(i))[in(i)]]cft + ^/QdL{t) 
ui(0) = 0, 



while u k , k = 2, . . . , n is the solution of 

\du k {t) = [Au k {t) + VF{<j){t))[u k {t)]dt + $ k (t)dt, 
\« fc (0) = 0. 

The remainder R n (t,s) is defined by 

n 

R n (t,e) := u(t) - 0(t) - J2 £k Mt) 

k=l 

ft / n ™ \ 

= / e^ A f («(s)) - - ^^Vf (^( S ))K( S )] - ^^^(s) ds, 

70 V fe=l fc=2 / 

and verifies the following inequality 



(ACP k ) 



(6.10) 



Sup |i?„(>,£)|^ 

te[o,r] 



<C p£ ' i+1 , 



wii/i a constant C p > 0. 



Proof. Let us define R n (t,e), n £ N, as stated in the theorem. Since by construction 

• 0(t) = e tA M° + /* e^-^ A F(<f)(s))ds (cf. Definition mj); 

• u(t) = e M u° + jle^- s ^ A F(u(s))ds + eL A {t) (cf. Definition SU); 

• Ul (t) = j t e i - t ~ s ') A \JF((j){s))[u 1 {s)]ds + L A (t) (cf. Proposition O and Definition 

• u k {t) = ^e^ A \7F(cj)(s))[u k (s)}ds + J* e^ A ^k(s)ds for k = 2,...,n, with $ fe (s) := 
$fc(</>(s)) [ui(s), . . . ,u fe _i(s)] defined in (|5.14|) (cf. Theorem IQl and Definition |4.6|) : 

we have 

ft f n n \ 

Rn(t, e) = / e^ A F(u(s)) - F^s)) - V s k VF^( S ))[u k (s)} - V £ fc <I> fc ( S ) ds . 

J ° V fe=l fc=2 / 
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Recalling that R[ n) (<j>(s); h(s), e) = F(u(s)) - F(^(s)) - £Li e k VF^(s))[u k (s)} - EL 2 £fc ^( s ) 
we get: 



sup \R n (t,e)\ p H 
te[o,r] 



<E 


sup 




te[o,T] 


< E 


sup 




te[o,T] 


< E 


sup 




te[o,T] 



;(*- s ) A i?["V(s);/i(s),e)* 



(t-s)A|ip 



?(")/ 



-w(t-«)P da 



(6.11) 



< p p(»+l) 

for some positive constant C njP (depending on n,p, but not on e), where in the second and third 
inequality we have used the contraction property of the semigroup generated by A. Now recalling 
Lemma I57T1 the inequality in Theorem 16.31 follows. 

□ 

Example 6.4. Our results apply in particular to stochastic PDEs describing the FitzHugh-Nagumo 
equation with a Levy noise perturbation (related to those studied with a Gaussian noise, for example, 
in [SBInJIZI] and [24]). 

The reference equation is given by (see [JH equation (1.1)]) 



(6.12) 



d t v(t,x) = d x (c(x)d x v(t,x)) -p(x)v(t,x) - w(t,x) + f(v(t,x)) +eL 1 (t,x), 
dtw(t, x) = jv(t, x) — aw{t, x) + eL2(t, x), 
d x v(t,0) = d x v(t,l) = 0, 
v(0,x) =v (x), w{0,x) = w Q (x), 

with the parameter e > in front of the noise, where u, w arc real valued random variables, a, 7 
are strictly positive real phcnomenological constants and c,p are strictly positive smooth functions 
on [0,1]. Moreover, the initial values vo,wq are in C([0, 1]). The nonlinear term is of the form 
f(v) = ~v(v — l)(t> — £), where £ € (0, 1). Finally L±,L2 are independent Qi-Levy processes with 
values in L 2 (0, 1), with Qi positive trace class commuting operators, commuting also with Aq, Aq 
being defined below. The above equation can be rewritten in the form of an infinite dimensional 
stochastic evolution equation on the space 

H :=L 2 (0, 1) x L 2 (0, 1) 

by introducing the following operators: 

A := d x c(x)d x , 

D(A ):= {ueH 2 (0, l);v x (0) = V x (l)} , acting in L 2 (0, 1) 

and 



(6.13) 



.4 



Ao-p -I 
7/ al 
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with domain D(A) := D(A ) x L 2 (0, 1), and 



F {l) = {~ V(V ) ' D{F) " L6(0 ' 1} X L2(0 ' 1} ' 

Further, we introduce the Banach space K := L ls (0, 1) x L 2 (0, 1), endowed with the norm | • \k := 
| • 1 6 + | • 1 2 and consider u° € K. In this way, the equation (|6.12[) can be rewritten as 

du(t) = Au{t) + F(u(t))dt + y/QdL{t) 
u(0) = u° := (v°,w°) e K 

with A and F satisfying Hypothesis 13.11 when £ 2 — £ + 1 < 3 min^gjo^i p{x). In fact, the properties 
of the two operators A and F can be determined starting from the problems considered in |24| and 
]26j . In particular from [551 Section 2.2] the estimates on the nonlinear term F and its derivatives 
can easily be deduced. Moreover we claim that the stochastic convolution 

L A {t) := f e^ A dL(s), 
Jo 

(where e tA 7 t > denotes the semigroup generated by A) is well-defined and admits a continuous 
version with values into the space K . This fact can be proved by an application of [63] and its proof, 
taking into account that the domain of fractional powers of A are contained in K (cf. Appendix A 
- in particular Example A. 5. 2 - in [32]) and moreover we are assuming TrQ < oo. 

Then by Theorem 16.31 we get an asymptotic expansion in powers of e > of the solution, in 
terms of solutions of the corresponding deterministic FitzHugh-Nagumo equation and the solution 
of a system of (explicit) linear (non homogeneous) stochastic equations. The expansion holds for all 
orders in e > 0. The remainders are estimated according to Theorem [63] These results should allow 
to obtain rigously results similarly to those obtained numerically up to second order in e in |69[ 172] 
in which the noise was of Gaussian type. Tuckwell, in particular, has made heuristic expansions up 
to second order in e for the mean and the variance of the solution process u = (u(i))t>o (see [69]l72]). 
proving in particular that one has enhancement (respectively reduction) of the mean according to 
whether the expansion is around which stable point of the stationary deterministic equation. 
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